Hyperbolicity of Hamiltonian formulations in 
General Relativity 



Ronny Richter 

Mathematisches Institut, Auf der Morgenstelle 10, 72076 Tubingen, Germany 
David Hilditch 

Theoretical Physics Institute, University of Jena, 07743 Jena, Germany 

Abstract. Well-posedness of the initial (boundary) value problem is an essential property, 
both of meaningful physical models and of numerical applications. To prove well-posedness 
of wave- type equations their level of hyperbolicity is an essential ingredient. We develop 
helpful tools and classify a large class of Hamiltonian versions of Einstein's equations with live 
gauge conditions with respect to their hyperbolicity. Finally we find a symmetric hyperbolic 
Hamiltonian formulation that allows for gauge conditions which are similar to the puncture 
gauge. 



1. Motivation 

Hamiltonian formulations play a crucial role in many areas of theoretical physics. Their key 
properties are the exact conservation of an energy and the symplecticity of the time evolution 
map. Certain properties of those systems can be obtained in numerical evolutions as well, using 
either symplectic [TJ or energy-preserving integrators [2J. 

Hamiltonian formulations for General Relativity (GR) were first considered in the 1950's 
OH]. Today the ADM formulation [5] is the most popular one. Yet, in numerical relativity, 
where one is interested in an initial value problem (IVP), Hamiltonian formulations are rarely 
used as the basis of evolution schemes. One can identify a reason for this in the mathematical 
properties of the ADM system and the numerical schemes used. 

The ADM system is composed of six evolution equations and four constraints. Furthermore 
it possesses a gauge symmetry. It can be shown that the constraint evolution system is closed, 
i.e. if the constraints are satisfied in the initial hypersurface then they are satisfied for all times. 

In numerical applications one usually prefers to solve for the constraints only in the initial 
data, because the alternative of solving the constraints in every timestep is computationally 
expensive. However, the IVP for the dynamical subsystem of the ADM equations is ill-posed. 
It is therefore not possible to build numerical schemes based on the dynamical ADM equations. 

A widely used method to avoid this problem and to deal with well-posed IVPs is to fix 
the gauge and to modify the ADM equations appropriately. The dynamical ADM equations 
are supplemented by evolution equations for the lapse and shift functions and multiples of the 
constraints are added to the dynamical ADM equations, such that the IVP for the full system is 
well-posed (the solution depends continuously on the given initial data). Here we are interested 
in the question how this technique can be applied while keeping the Hamiltonian structure. 



2. Hyperbolicity 

In order to construct systems with well-posed IVPs one relies on notions of hyperbolicity. It can 
be shown that strongly (symmetric) hyperbolic systems possess well-posed initial (boundary) 
value problems (of course one needs appropriate boundary conditions in the case of IBVPs) 
[61 [7]. The advantage of this approach is that the proof of hyperbolicity involves essentially 
linear algebra. Strong and symmetric hyperbolicity are defined as follows [8]. 
A first order in time, second order in space system of the form 

d t v = A[diV + Aiv + A 2 w + a, (la) 
d t w = B\ j did jV + B\d iV + B x v + B\d iW + B 2 w + b. (lb) 

is called strongly hyperbolic if the principal symbol 

ps _( A\ Si A 2 \ 
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has a complete set of eigenvectors (with real eigenvalues) that depend continuously on the spatial 
vector Si. 

The system ([1]) is called symmetric hyperbolic if there exists a conserved positive definite 
energy, E = j edx [8] of the form 

We denote H % i a symmetrizer. Furthermore if there is a matrix H % i such that u\H^Uj 
is a conserved quantity, but not necessarily positive definite then we call a candidate 
symmetrizer. 

3. Hamiltonian formulations 

A Hamiltonian formulation is given by the specification of a Hamiltonian H(q,p) for the 
system, constructed from canonical positions and momenta (q;p). Here the Hamiltonian can be 
expressed as the integral over space of a Hamiltonian density 

U = / d 3 xU D (q,diq,p,dip), (4) 

where the Hamiltonian density is a local function of the canonical variables and their spatial 
derivatives. One obtains the following canonical equations of motion 

_ &K D dn D , e , 

dU D dn D v 

dtP= -^q- +d ^r (5b) 

A system that can be written in the form ([5]) has Hamiltonian structure. 

In order to keep the Hamiltonian structure while modifying the evolution system a good 
strategy is to modify the Hamiltonian density appropriately and derive the equations of motion 
afterwards. The key ideas in this approach were presented in [9]. 

One starts from the ADM Hamiltonian density 



^adm = / {-aC + 2pC l )^d 3 x, , 



(6) 



where a is the lapse function, (3 l is the shift, 7 is the determinant of the 3- metric and C, Ci are 
scalar and vector constraints respectively. The constraints are 



C = - ~7rVj - hR, Ci = -2%Ay' fc , (7) 

with the 3-metric "fij and its canonically conjugate momentum ir l;) . The dynamical degrees of 
freedom are (7^ , n % i ) . 

Then in order to provide evolution equations for lapse and shift one needs to include them 
in the set of dynamical degrees of freedom. This also means that they need conjugate variables, 
which we denote a and pi respectively. One finds that these momenta vanish for solutions of 
GR, i.e. <7 = and pi = 0. 

After modifying the phase space we must also work with a different Hamiltonian. One adds 
to the ADM Hamiltonian terms that allow for appropriate gauge conditions, but do not change 
the physics. We use 

n = Uadm + Ughg + J d 3 x (Act + Wfa) , (8) 

where Hadm + Hghg is the Hamiltonian for Brown's generalized harmonic formulation [9]. The 
terms that appear in the principal part are 

Ughg = /3VA« + Fpjdifi + a 2 Pl T) k ^ k - apWia 

+ ( - 4a 2 7i j vr i V + aV - lo'V/',- • ( 9 ) 
°V7 

The terms in A and that affect the principal part of the equations of motion are linear in the 
canonical momenta and in the first spatial derivatives of the positions. With the restrictions 
described in |10| we come to the following form 

A = -C ia 2 7" 1/2 7y^ + C A aS~ 1/2 <? + C 7 (aDifi 1 - ^a^^d^A (10) 

and 

W = C 2 a 2 T) kl jk + C 3 a 2 T k kjl ji - C 5 af j D ja - C^~ l/2 a 3 Pj Y j . (11) 

With this choice the principal symbol depends on the shift in a trivial way, which simplifies the 
hyperbolicity analysis. 

The consequence of using this approach is that one cannot choose equations of motion for 
lapse and shift and constraint additions independently. In fact, if the gauge conditions are fixed 
then constraint additions are fixed as well. 



4. Crucial Techniques in the hyperbolicity analysis 

Given the Hamiltonian ([8]) with parameters Ci it is straightforward to compute the corresponding 
equations of motion and their principal symbol ([21). 

As we describe in [lOj . since the principal symbol depends only trivially on the shift, it has 
essentially the following structure 

p- = (° Y x )- m 



This matrix decomposes further into scalar, vector and trace-free tensor blocks which simplifies 
the analysis further. The most complicated part in the proof is then the scalar block. The 
corresponding sub blocks of X and Y both have upper block triangular form: 

with 2x2 matrices A, B and C. This structure can be used to reduce the analysis completely 
to 2 x 2 matrices. Finally one can identify all three families of strongly hyperbolic formulations 

Concerning symmetric hyperbolicity we consider only the strongly hyperbolic formulations, 
because symmetric implies strong hyperbolicity. The first step is to construct candidate 
symmetrizers. This can be achieved by making an ansatz for the candidate and solving linear 
equations for the parameters in that ansatz. 

Having a set of candidate symmetrizers one can ask for positive definite matrices in that set. 
This is the most complicated step in the analysis and we are not able to completely classify the 
formulations with respect to their symmetric hyperbolicity. 

We can however identify strongly hyperbolic formulations that are not symmetric hyperbolic 
and have also found 2-parameter families of symmetric hyperbolic formulations. 

Thus far we have treated the formulation parameters Cj as constants. But one can view them 
as functions of ("fij, a, /3 % ) as well, which effects neither the physics nor the level of hyperbolicity, 
but does change the equations of motion. One can use this freedom to derive a formulation that 
includes a gauge condition very similar to the popular puncture gauge. The equations of motion 
for lapse and shift are 

d t a = (3 l dia - ixlc?K + —^—fi L a 3 a, (14) 

dt? = pdjP + M s7 1/3 r} fc 7 ife + ^(M57 1/3 " o?)Y\^ - vsj- 1/e ap l - aD*a. 

which is (near the puncture (a —> 0) and if the constraints a = and Pi = are satisfied) very 
close to the puncture gauge condition 

d t a = p i d i a- f i L a 2 K, = fidft + /x 5 F - r,/?, (15) 

where 

r = 7 1/3 r^7 ife + ^V'rV (16) 

Acknowledgments 

This work was supported by DFG grant SFB/Transregio 7 "Gravitational Wave Astronomy". 
[1] E. Hairer, C. Lubich, and G. Wanner. Springer Series in Computational Mathematics 31. Berlin: Springer., 
2006. 

[2] E. Hairer. J. Numer. Anal. Ind. Appl. Math. 5 (2010) 73-84 
[3] F.A.E. Pirani, A. Schild, R. Skinner. Phys. Rev. 87 (1952) 452 
[4] P. Dirac. Proc. Roy. Soc. Lond., A 246:333, 1958. 

[5] R. Arnowitt, S. Deser, C.W. Misner. in L. Witten (Ed.), Gravitation: An Introduction to Current Research. 
Wiley, New York (1962) 227 [arXiv:gr-qc/0405109] 

[6] B. Gustaffson, H. Kreiss, J. Oliger. Wiley, New York 1995 

[7] G. Nagy, O.E. Ortiz O.A. Reula. Phys. Rev. D70 (2004) 044012 [arXiv:gr-qc/0402123] 

[8] C. Gundlach, J.M. Martm-Garcfa. Class. Quant. Grav. 23 (2006) S387-S404 [arXiv:gr-qc/0506037] 
[9] J.D. Brown arXiv:gr-qc/0803.0334 (2008) 
[10] D. Hilditch, R. Richter. arXiv:gr-qc/1002.4119 (2010) 



